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Abstract 

A thermodynamically consistent model of non-classical coupled non-linear ther- 
moelasticity capable of accounting for thermal wave propagation is proposed. 
The heat flux is assumed to consist of both additive energetic and dissipative 
components. Constitutive relations for the stress, the entropy and the energetic 
component of the heat flux are derived in a thermodynamically consistent man¬ 
ner. A Lyapunov function for the dynamics is obtained for the case in which 
the surface of the continuum body is maintained at a reference temperature. It 
is shown that the system is non-linearly stable. The linearized model is shown 
to be similar to the type III model of Green and Naghdi, except for some minor 
differences in the interpretations of some of the parameters. 
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1. Introduction 

The propagation of heat energy in a rigid body is governed by the energy 
balance equation which reads 


e = —divq + r. (1) 

Here e, q, and r denote the internal energy, heat flux, and heat source respec¬ 
tively. The superposed dot denotes time derivative and div[ • ] divergence. In the 
classical theory of heat conduction with Fourier’s law [lj , the balance equation 
m is supplemented by a constitutive equation for the heat flux q satisfying the 
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inequality 


q • V0 < 0, 


( 2 ) 


where 0 is the absolute temperature. The inequality ©, also referred to as heat 
conduction inequality , states that heat only flows by the mechanism of diffusion 
from regions of higher to those of lower temperature distributions. In other 
words, the vector quantity describing the heat flow, that is, the heat flux q, 
always points in the direction of decreasing temperature distribution. 

However, for example, in the linear case this mechanism leads to a parabolic 
partial differential equation which is characterised by infinite speed of prop¬ 
agation of localized thermal disturbances, a paradoxical phenomenon from a 
physical point of view. 

Despite this non-physical prediction, the classical theory of heat conduction 
by Fourier’s law has been successful for a broad range of engineering applica¬ 
tions. However, as length scales decrease, this phenomenon of instantaneous 
propagation of thermal disturbances becomes more dominant at low ranges of 
temperature near absolute zero, so that approximations used in the classical 
theory of heat conduction lose validity Q • 

The earliest known conjecture on the existence of thermal propagation as 
waves, also known as the second sound phenomenon, was given by Nernst Q in 
1917. Later, in 1938 Tisze [ij and Landau [1] in 1941 independently suggested 
the possibility of thermal waves in superfluid liquid helium, at temperatures 
below the so-called lambda transition near 2.2 K [§j. Peshkov Q reported 
the first experimental evidence for the existence of second sound in 2He. In 
his work, Peshkov suggests that second sound might also be observed in pure 
crystalline materials on the basis of similarities between crystal materials and 
liquid helium. Laser pulsing experiments have shown that second sound can 
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propagate in high-purity crystals of 4He @], 3He [9], NaFfloj , and Bi 

Cattaneo [12| was the first to introduce a non-Fourier theory of heat con¬ 
duction in order to overcome the paradoxical prediction of the classical Fourier 
theory. His work is based on the concept of relaxing the heat flux from the 
classical Fourier law to obtain a constitutive relation with a non-Fourier effect. 
There have been several other attempts to develop continuum theories capable 
of predicting thermal waves propagating at finite speeds for various types of me¬ 
dia. Among these, the works in 13-17] are notable. The constitutive equation 
for the heat flux according to Cattaneo is examined in terms of thermodynamics 
in nnm. Later this result was extended in [20j to the case where deformation 
is applicable. 

A relatively more recent theory of non-classical heat conduction with and 
without deformation was proposed by Green and Naghdi 21 24]. Their work is 
based on the introduction of three types of constitutive relation for the heat flux, 
thereby resulting in three different models, namely type I, which is the classical 
theory, type II, a purely hyperbolic model which allows the propagation of a 
heat pulse without damping, and type III, which is the combination of the first 
two. 

In recent years there has been a considerable amount of interest concerning 
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the theory of Green and Naghdi. 
be found in 


25^271 . 


An extensive overview of the theory can 
Theoretical results addressing existence and uniqueness 
and exponential stability (30j | have also been investigated for some types 
of the theory. The designing of appropriate numerical methods has also been 
addressed in 


3l|, |32|. 


This work is concerned with a thermodynamically consistent formulation 
of a fully non-linear coupled problem of non-classical thermoelasticity inspired 
by that of Green and Naghdi. The formulation is based on the basic laws 
of continuum thermodynamics, the balance laws of momentum, the balance 
of energy, and the entropy imbalance. However, the point of departure from 
the classical theory comes from two assumptions: the first is that the heat 
flux is additively composed of two parts, namely the dissipative and energetic 
components, and the second is that a material derivative of a time primitive of 
the absolute temperature is assumed to be the proportionality constant of the 
heat and the entropy conjugate pairs. 

Thermodynamic restrictions on the constitutive relations are derived using 
the procedure of Coleman and Noll 33]. Stability of the system of partial differ¬ 
ential equations governing the thermomechanical coupling in the non-classical 
regime is proved in the sense of Lyapunov. The other notable aspect of the 
model is that the linearized theory is similar to that of by Green and Naghdi 
except for some differences in the interpretation of the material parameters. 

The rest of the paper is organised as follows. In Section [2] geometric and 
kinematical descriptions of the continuum body are presented. In Section [3] the 
non-classical model describing the coupling of mechanical deformation and non- 
classical heat conduction is formulated based on the laws of thermodynamics. 
Constitutive relations for the stress, the entropy and the energetic component 
of the heat flux are derived from a free energy via the Coleman-Noll procedure. 
Next, in Section [4] the initial boundary value problem (IBVP) of non-classical 
thermoelasticity is summarised. A class of physically meaningful initial and 
boundary conditions are also proposed. A Lyapunov function for the dynamics 
generated by the IBVP rendering the system non-lincarly stable is obtained in 
Section [5] Finally in Section [6j the linearised form of the IBVP is summarised. 


2. Kinematics: Lagrangian description 

Consider a continuum body occupying the reference configuration fl C R d , 
where d = 1, 2, or 3, with smooth boundary T and material point denoted 
by X G Cl. The surface outward unit normal in the reference configuration is 
denoted by N. The motion of the body is parametrized by a C 1 orientation¬ 
preserving map: 

<p(X,t) :I->R d , X € IT, (3) 

where fl represents the closure of fl, or the union of the interior of f l and its 
boundary T, and I = [0,T] denotes the time domain. The displacement and 
velocity fields in the Lagrangian description associated with a particle located 


3 









at point X at time t € I respectively are 

u(X,i) = ¥>(X,i)-X, v(X,i) = ^u(X,i). (4) 

In what follows, partial time derivatives will be denoted by superimposed dots 
so that b = do/dt for some field o in the material description. The deformation 
gradient F, the primary measure of deformation, is defined by 

F = V« + 1, (5) 

where the operator V denotes the gradient with respect to the reference con¬ 
figuration, and 1 is the identity second-order tensor. For an arbitrary material 
point X and time t S I, an infinitesimal reference volume element dfl associated 
with the material point is mapped into a volume element d u in the current 
configuration by the Jacobian J(X,t) = detF(X,f) via the relation 

duj = J{X,t)dfl. (6) 

Thus, an admissible configuration is one in which the Jacobian J is positive. 

3. Non-linear coupled thermo-mechanics in the non-classical frame¬ 
work 

3.1. Balance laws 

• The local forms of the momentum balance are 

p 0 ii = DivP + p 0 b , (7) 

PF t = FP t . (8) 

Here po , b and P are, respectively, the density in the reference configura¬ 
tion, the body force density and the first Piola-Kirchhoff stress tensor. The 
superscript T denotes the transpose operation on second-order tensors. 

• The balance of energy (first law of thermodynamics) is given by 

£ = P : F - DivQ + R, (9) 

where £, Q, and R are the internal energy, heat flux, and heat source. 

• The entropy imbalance (second law of thermodynamics) is given by the 
inequality 

i) > —DivH + S, (10) 

where r/, H, and S are the entropy, entropy flux vector, and entropy source, 
respectively. 
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The absolute temperature 0 is used to relate heat flux (resp. heat source) with 
entropy flux (resp. entropy source) in the classical theory (see, for example [34 j 
pp. 187). In our formulation, motivated by the theory of Green and Naghdi, 
we additionally assume that the existence of a scalar field a, referred to as the 
thermal displacement , such that 


and that 


Q 


a 


a > 0, 

(11) 

and S = —. 

a 

(12) 


We then proceed to define the absolute temperature by 


0 = a. 


Next, substitution of m into (in leads to 
fj > —Div 


^ or 0?? > -DivQ + • V0 + R. (13) 


By subtracting the inequality (fl3l) from the energy balance equation © and 
subtracting Qr/ from both sides of the resulting inequality, we obtain 


^[£-0??] <P :F- 07 7 -iQ V0. 


(14) 


On defining the Helmholtz free-energy through the Legendre transformation 


= £ - 077, (15) 

(ED yields the local form of the free-energy imbalance 

4' - P : F + 0r ? + • V0 = -05 < 0. (16) 

Here 3 denotes the rate of entropy production (while 03 is the rate of energy 
dissipation). 


3.2. Constitutive theory 

As mentioned above in the motivation, we assume that the heat flux Q is 
split additively as 

Q = Q e + Q d- (17) 

We refer to and Qd as the energetic and dissipative components of the heat 
flux Q. As a consequence of this split, the free-energy imbalance (fill becomes 

'L - P : F + 07? + ^Q e • V0 + 07? + • V0 = -03 < 0. (18) 
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Looking at the variables involved in the free-energy imbalance m , we as¬ 
sume that the state of the system under consideration depends on the set of 
state variables 

si = {F,0, Va, V0}. (19) 

Thus the free energy the Piola stress P, the entropy 77 , and the heat flux 
components Q e and Qd are determined by constitutive equations of the form 


T = P = P (si), rj = i](si), 

Q E = Q E(si), Q D = QdM- (20) 


We now apply the Coleman-Noll procedure to ensure that the constitutive equa¬ 
tions (l20l) satisfy the laws of thermodynamics. If the constitutive equation for 
the free energy is differentiated with respect to time, one obtains 


• ■ 3tf ■ d<S> ■ 3$ ■ 

Vp = - • F -I-0 -I-- A H-- G 

3F ' 30 3A 3G ’ 


( 21 ) 


where A = Va and G = V0. 

Substitution of the time derivative of the free energy m into the free im¬ 
balance m yields 


d'L 

3F 


-P :F + 


3T 

30 


+ 77 0 + 


3'L 

3A 


1 

0 


Q E A 


3^ 

3G 


■G + j-Qd • V0 < 0, 


0 


( 22 ) 


which must be satisfied for all states. 

Since F, 0, A, G, and their time primitives can be chosen arbitrarily in 
order to maintain the inequality (1221) we may choose the constitutive equations 
and thermodynamic restrictions as 


P 


3T _ 3T 1 
3F ’ ' n ~-~dQ 1 0 Qe 


34- 

3A’ 



and Qo-V0<O. (23) 


As consequence, the rate of entropy production in any thermodynamically ad¬ 
missible process becomes 


= - 02 ^- Ve - 


Note that A = V0. Equation (E51) 1 reveals that the free energy does not depend 
on V0. 

For instance, a class of such free energy functions in the non-classical case 
has the form 

’£ = T c + -KiVa • Va, (24) 

where T c denotes any classical free energy function for thermo-hyperelasticity 
and Ki a symmetric and positive-definite non-classical conductivity second- 
order tensor. Equation (l24l) and the constitutive equation (l23lE imply that the 
energetic component of heat flux Q# is linear in Va. 
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3.3. Field equation for the thermal conduction 

Consequently, by ED and ED. the time derivative of the free energy be¬ 
comes 

^ = P:F-r,0-^Q E -A. (25) 

Next, by the Legendre transformation ED, from (l25l) the time derivative of 
the internal energy is 

£ = 07) + P:F^iQ B .A. (26) 

By combining © and (l26l) we obtain the field equation of thermal conduction 
in the conservation form 

0t)= -DivQ+^Qs- VO + R. (27) 

From ED i ,2 an d assuming sufficiently smoothness on 4' with respect its argu¬ 
ments, 

dij d 2 ^ d 2 ^ _ <9P 

5F - ~ OFdQ ~ dQdF ~ 30' ( ’ 

Using equation ED the conservation form (1271) yields the field equation 

1 9P • 

p o c F 0 = —DivQ + —Q E ■ V0 + 0— : F + p 0 R , (29) 

where c F is the heat capacity and is defined by 


Poc F 


= -0 


d 2 ^ 

a© 2 ' 


(30) 


4. The initial-boundary value problem 

In this section we summarise the initial boundary-value problem for coupled 
non-classical thermoelasticity in a general non-linear framework. 

4-1. Local balance laws and constitutive equations 

The displacement-velocity relation, the balance of linear momentum J7D , the 
thermal displacement-temperature relation along with the balance of energy 
in the conservation form (1271) lead to the first-order (in time) system of non¬ 
linear partial differential equation governing the strong coupling of thermoelastic 
behaviour in the non-classical regime, given by 


u = v 

p 0 v = DivP + p 0 B 
d = 0 

0 ?) = —DivQ + jtQe • V0 + R 

> 


in llxl. 


(31) 
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The constitutive equations are derived from a potential denoting the classical 
Helmholtz free energy function, such that 


<94- <94> 

dF’ V ~~dQ 


Q = Q e + Q 


«9T 

d and Q e = — 

a A 


(32) 


Equation (|32l) together with a constitutive relation for the dissipative component 
Qd of the heat flux satisfying an inequality of the form Qd • V0 < 0 complete 
the set of constitutive equations for system m- 


4-2. Initial and boundary conditions 

Physically meaningful initial and boundary conditions in the non-classical 
theory are proposed in 32| . To define initial conditions, we assume an initial 
configuration in which the thermal displacement field is homogeneous. Thus, 
we use initial conditions of the form 


= v 


^li=0 


and 0| t=o = 0° 


(33) 


where (•)° are initially prescribed functions on Q. 

For boundary conditions, we first consider the decomposition of the bound¬ 
ary T into two mutually disjoint partition sets (T u , T t }, and (T e , T Q } such 
that 

r„ n r T = r e n r Q = 0, and r„ u r T = r e u r Q r. 

Hence physically meaningful boundary conditions have the form 


u = u 

on P 

u 

x I, 

and PN 

= T 

on l’ T 

x I, 

e = e 

on r e 

X I, 

and Q • N 

= Q 

on r Q 

X I, 


(34) 


where (•) denote prescribed functions on the corresponding partitions of T. 

Therefore, the system m together with the constitutive relations (1321) . the 
initial and boundary conditions (1531) and (1M1) conclude the local description of 
the initial-boundary value problem (IB VP) for coupled non-classical thermoe¬ 
lasticity. 


5. Stability of the continuous problem. The Lyapunov function 

Assume that there are no heat sources (R = 0) and the part r e of the 
boundary is fixed at a constant reference temperature 0 O the temperature of 
the medium in which the motion of the solid takes place. Moreover, a thermally 
homogeneous boundary condition (Q • N = Q = 0) is prescribed on T Q , and the 
body and boundary mechanical loadings are conservative; that is, denoting the 
potential energy of the external loadings by n ext , we have 

^n ext = - J p 0 B ■ v dn - j t • t, dr. (35) 
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Defining the relative temperature by d = 0 —@o, the functional for the dynamics 


E = [T + dr/ + • v] dfl + Ilext, (36) 

corresponding to the specific conditions just described above defines a Lyapunov 
function for the IB VP (IdlT) . This can be shown by differentiating E with respect 
to time as 


-J-E = 

d t 


n 


P : F - 77 © - — Q e • V0 + d ?7 + d?) + pqV ■ v 


dfl + -j-n ex t 
at 


0, 


= J ■ V 0 <m 

< 0 . 


(37) 


In the calculation leading to (1371) we make use of the following: the time deriva¬ 
tive of the free energy (l?5l) . the linear momentum balance 0 and (1551) together 
with integration by parts, and then the energy balance in conservation form 
(1271) with integration by parts, and finally flux partitions and the constitutive 
restriction (1231k on the dissipative component Q d of the heat flux. 

Equation J37J implies that the positive definite functional E is monotonically 
decreasing along the flow defined by the IBVP, which proves the non-linear 
stability of the continuous problem for the specific case outlined above. A 
similar result for the case of more general boundary conditions is not known, 
even in the classical case (see also 1351] h 


6. The linearized theory 

The IBVP (1311) is linearized about the reference configuration for which u = 
0, v = 0, a = 0 , and 0 = 0 o are assumed to be natural for the configuration. 
As a consequence, the linearized version of IBVP (l3lT) is given by 


u = v, 

pv = diver + p 6 , 

a = 0 , 

0 O 77 = — divq + r, 


(38) 


with constitutive equations 


di/j dip dip 

<T= Se’ r? = _ a0’ q = qe + qd ’ ge = -0o M’ and Qd 


where ip is the quadratic free energy function defined by 


ip = -£ : Ce — dm : e -77-d 2 H — KiVa • Va. 

2 20 o 2 


—K 2 V0, 

(39) 

(40) 
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Then, by (1391) and (HOl) . we obtain 

er = Ce — dm, 77 =^-'# + m : e, and q = — 0oKiVa — K 2 V0. (41) 

0o 

Here, <x, q 1 e, m are the stress tensor, the heat flux, the infinitesimal (sym¬ 
metric) strain, and thermomechanical coupling tensor, respectively. Under the 
assumption of isotropy, the fourth-order elasticity tensor C, and the non-classical 
conductivity tensor Ki are both symmetric and positive-definite, while the clas¬ 
sical heat conduction tensor K 2 is symmetric and positive-semidefinite. 

Note the similarity between this linearized model and the type III model 
by Green and Naghdi (see for example [32j]). The only difference is that the 
non-classical conductivity tensor Ki is replaced by @oKi. However, the model 
due to Green and Naghdi is not consistent with the laws of thermodynamics in 
the non-linear case. 

7. Conclusion 

A thermodynamically consistent formulation for coupled thermoelasticity at 
finite strain in a generalized framework was developed. The model encompasses 
the classical model based on Fourier’s law of heat conduction as well as the non- 
classical one in which a non-Fourier type heat conduction equation allows for the 
propagation of thermal energy as a wave (with finite speed). The formulation 
is based on two key postulates: the first is that the introduction of a scalar 
internal state variable, the thermal displacement, defined as the time primitive 
of the absolute temperature, and the second is that the heat flux is composed 
of both energetic and dissipative components. Physically meaningful initial 
and boundary conditions were presented. A Lyapunov type functional for the 
governing IBVP was suggested for some class of boundary conditions. It was 
shown that the system is non-linearly stable. The linearized model was shown 
to be similar to the type III model of Green and Naghdi. 

The design and analysis of efficient numerical algorithms for the non-linear 
thermomechanical model will be described in future work. For such non-linear 
and strongly coupled problem, it seems natural to use of shock capturing schemes 
such as the space-time discontinuous Galerkin approach together with an operator¬ 
splitting technique to effectively decouple the system of governing equations 
without violating the stability of the full system. 
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